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We investigate the connection between the geometry of Favoured Local Structures (FLS) in liquids 
and the associated liquid and solid properties. We introduce a lattice spin model - the FLS model 
on a face-centered cubic lattice - where this geometry can be arbitrarily chosen among a discrete 
set of 115 possible FLS. We find crystalline groundstates for all choices of a single FLS. Sampling 
all possible FLS’s, we identify the following trends: i) low symmetry FLS’s produce larger crystal 
unit cells but not necessarily higher energy groundstates, ii) chiral FLS’s exhibit to peculiarly poor 
packing properties, iii) accumulation of FLS’s in supercooled liquids is linked to large crystal unit 
cells, and iv) low symmetry FLS’s tend to find metastable structures on cooling. 


I. INTRODUCTION 

For a given interaction potential between particles, 
some local arrangements will be particularly stable and 
these favoured local structures (FLS) will thus accumu¬ 
late in the liquid upon cooling. There is a considerable 
body of literature on the geometric and kinetic conse¬ 
quences of these FLS’s in supercooled liquids, which have 
been the subject of a number of recent reviews mm- A 
possibility that has remained of perennial interest since 
first suggested by Frank in 1952 [3] is that if a single 
FLS dominates the liquid structure and its geometry is 
incompatible with crystalline symmetries, it could sta¬ 
bilize the liquid down to the point of dynamical arrest. 
This is the starting point of the geometrical frustration 
approach to understanding the glass transition [4]. The 
fundamental question of this field can be posed as fol¬ 
lows: might there be favoured local structures of such 
geometrical awkwardness that they either do not have a 
periodic ground state or, if they do, it is so high in en¬ 
ergy that the freezing point is forced down well below the 
limit of kinetic accessibility? 

To properly address this question it is important that 
we consider the consequences of any particular set of 
FLS’s for both the crystal phase as well as the liquid 
state. Indeed, the stability of the liquid depends cru¬ 
cially on the stability of the crystal and the kinetics of 
transformation into that ordered structure. In adopting 
this perspective, physical consistency requires that the 
crystal structure must be assembled out of the same set 
of FLS’s available to the liquid. The fact that the sim¬ 
ulation studies of liquid FLS’s typically report little or 
no concentration of the crystalline local structure mm 
(despite their stability) is, we suggest, a consequence of 
the high symmetry of the local environments found in the 
small unit cell crystals frequently represented in simula¬ 
tion models. As we have previously shown [5], a liquid 
can only accommodate a low concentration of high sym¬ 
metry FLS’s due to their high entropy cost. In cases 
where the crystal unit cell is larger than just the nearest 


neighbour coordination, the crystal local structure has 
been observed at high densities in the supercooled liq¬ 
uid [6]. 

Reports of systematic searches of low energy or high 
density crystal structures have appeared with increas¬ 
ing frequency in recent years m- These studies have 
highlighted the extraordinary variety of crystal structures 
that can be accessed by the parametrization of simple 
shapes. The role of local structure in liquids remains a 
more challenging question. The structure of molten salts 
has generated an extensive literature [12] as has the struc¬ 
tures in liquids characterized by local tetrahedral order 
such as silica m- Even in cases when the local structures 
are complex and multiple, there exists a strong conviction 
that they represent an essential key to understanding the 
properties of the low temperature phases. A considerable 
body of work addresses the favoured local structures in 
binary and ternary atomic alloys [2]. These local struc¬ 
tures have been associated with dynamic arrest and slow 
crystallization rates [2|. 

Our objective is to complete a systematic study of how 
changes in the geometry of the FLS influence the prop¬ 
erties of the liquid and crystal phases. To carry this out 
through adjustments of the interparticle potential would 
be difficult, requiring the invention of a set of complex 
potentials for the task. In this article we shall circum¬ 
vent this problem by explicitly defining our Hamiltonian 
in terms of the FLS’s themselves. In this way we simply 
assign potential energies to each of the possible FLS’s 
so as to favour selected local structures. As the number 
of possible FLS’s can be large, even this simple recipe 
can generate an enormous number of possible Hamiltoni¬ 
ans so we shall restrict our attention to the cases where 
only a single local structure at a time is favoured. With 
this restriction we shall not be able to directly address 
the original Frank picture mm which requires at least 
two FLS’s - an icosahedral structure assumed to domi¬ 
nate the liquid structure and the cube octahedron that 
characterises the local structure of the face centred cubic 
(fee) crystal. The influence of multiple FLS’s we leave 
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for future work. 

In this paper we implement this FLS Hamiltonian us¬ 
ing Ising spins on the vertices of an fee lattice. We have 
previously reported on the features of this FLS model 
in two dimensions, where we explored the relationship 
between the FLS and the crystal structure [7] , liquid en¬ 
tropy [5] , the freezing transition [8] and the role of chiral- 
ity [9] . In this paper we shall examine how these various 
observations are influenced by the increase to 3D. The 
immediate challenge of the 3D FLS model is the signifi¬ 
cant increase in the number of possible FLS’s - from 9 in 
2D to 115 in 3D. 

The paper is organised as follows. In Section [TTj we 
detail the definition of the three-dimensional model. We 
then turn to the description of the different regimes of 
interest: in Section m we investigate the structural 
properties of the ground states, finding that they are all 
crystalline but can be very complex. In Section |IV| we 
describe the thermal behaviour of the FLS model: liq¬ 
uid phase, supercooling and freezing transitions. Finally, 
Section [V| is devoted to an analytic study of the pre¬ 
viously studied quantities, using statistical methods to 
connect structural and thermal properties. We discuss 
the consequences of our findings in Section [VH Details of 
the groundstate properties and transition behaviour for 
each FLS is provided as Supplementary Material HO]. 


II. THE MODEL 

The Favoured Local Structures model is a model of 
Ising spins (ie. binary variables) located at the vertices 
of a regular lattice. In this article we focus on the case 
of a three-dimensional face-centered cubic (fee) lattice. 
The local structure at a given lattice site is defined as 
the spin conformation of the nearest neighbours of this 
site (regardless of the value of the spin at said site). Thus 
there are 2^ possible local structures at a given site, where 
2 : = 12 is the connectivity of the lattice. Among this large 
number of structures, some differ only by a rotation and 
are treated as indistinguishable, so that this model is 
isotropic. An instance of the model is obtained by select¬ 
ing one local structure and making it favoured by directly 
stabilizing it. In order to do so, we attribute an energy of 
— 1 to sites whose local environment corresponds to this 
favoured loeal strueture (up to a rotation), while all other 
sites have zero energy. The total energy of a spin config¬ 
uration is thus simply, up to a minus sign, the number of 
sites lying in the FLS. We consider configurations of such 
systems on a large lattice with periodic boundary condi¬ 
tions, in the canonical ensemble. The temperature thus 
controls the density of Favoured Local Structures, and 
the zero-temperature limit corresponds to dense packing 
of FLS on the lattice. 

On the two-dimensional triangular lattice previously 
studied EllZHS], z = 6 and thus there were 64 possible lo¬ 
cal structures, falling into 14 rotationally distinct classes. 
As the model turns out to be the same when one FLS or 


its spin inverted variant is selected, there were 9 possi¬ 
ble choices left. Finally, two of these local structure were 
ehiral^ i. e. they lacked a reflection symmetry, and formed 
an enantiomer pair. While the model would be identical 
by selecting either of those enantiomers, we found inter¬ 
esting to distinguish the case when only a one of them 
is favoured - the ehiral case - from that when both are 
favoured - the raeemie case. Physically, the racemic case 
might correspond either to an equimolar mixture of enan¬ 
tiomers or to a flexible molecule with two enantiomeric 
conformations. This left us with 9 possible choices for 
the FLS, all of which were studied in detail. 

In this article we consider the case of a face-centered 
cubic lattice, for which each site has 2 : = 12 neighbours. 
The 4096 local structures reduce to 115 after removing 
those which are related by rotation or spin inversion of 
another. Among them, there are 39 enantiomer pairs for 
which we distinguish the racemic and the chiral cases, 
while the remaining 37 structures are achiral. A selec¬ 
tion of these structures, where spins are represented by 
coloured spheres, is presented in the insets of Figure 
All these structures are depicted, labelled and individu¬ 
ally studied in the Supplementary Material m- We now 
report on the study of these 115 possible systems. 

The main question we want to address is the influence 
of the local geometry on the macroscopic properties of 
the system. To quantify the degree of symmetry of the 
local structure, we introduce the size S of the spatial 
point group of an FLS. This number ranges from S = 1 
for those FLS’s with no symmetry save the identity op¬ 
eration to 5 = 48, the size of the complete symmetry 
group of the lattice. The latter is attained only by one 
FLS, the all-up environment labelled 0 , corresponding 
to ferromagnetic ordering. The inclusion of the racemic 
mixtures complicates the test for the influence of 5, since 
the added flexibility due to favouring two distinct enan¬ 
tiomers is not reflected by this symmetry parameter. Be¬ 
sides, we show in Section [VB] that the chiral symmetry is 
peculiar in our model. For these reason, we distinguish 
the sets of achiral, chiral and racemic FLS throughout 
the paper. 


III. STRUCTURAL PROPERTIES AND 
GROUND STATES 

In this Section, we present the ground state properties 
of the model, corresponding to dense packing of the FLS 
on the lattice. 


A. Ground state energies 

The ground state energies corresponding to minus 
the maximum packing density of the FLS, have been de¬ 
termined for all choices of FLS. These energies and the 
associated crystalline structures were determined by an 
enumerative process similar to that introduced by Hart 
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FIG. 1: A selection of ground state structures for various choices of the favoured local structure. This 
representation shows a portion of an infinite crystal, cut along different planes. The associated FLS is depicted in 
inset, while the labels correspond to those employed in Supplementary Material. For chiral FLS’s, the letters c and 
r respectively correspond to systems with one and both enantiomers favoured. 


and Forcade m, detailed in the Appendix. We found no 
evidence of a non-crystalline structure more stable than 
the crystals identified by this method. A selection of crys¬ 
talline structures is presented in Figure A complete 
presentation of these structures is provided in ref. uni. 


The distribution of ground state energies is plotted 
in Figure These energies are integer fractions, cor¬ 
responding to minus the ratio of the number of FLS in 
the crystal cell to the size Z of this crystal. The ground 
state energies range from —1 (corresponding to all sites in 
the FLS, attained by three local structures) to —1/4, the 
lowest packing fraction obtained (attained by five chiral 
FLS’s). Notably, chiral FLS’s tend to have higher ground 
state energy than other structures, a fact that will be dis¬ 
cussed in Section |VB| There is a particularly large peak 
at Eo = —1/2 characterized mostly by structures consist 
of an alternating 2D layers that are either completely 
filled or completely empty of FLS. 


With an average of (Eo) = —0.51(0.17) [30], the mini¬ 
mal energies are generally considerably higher than those 
observed in 2D [ 7 ], for which (Eo) = —0.73(0.13). This 
effect of dimension is analogous to the well documented 
decrease in maximum packing density for hard spheres 
with increasing dimensionality m- The increase in the 
number of nearest neighbour sites resulting from an in¬ 
crease in spatial dimension appears to generate more con¬ 
straints with regards to packing than it provides degrees 
of freedom to resolve them, giving rise to the observed 



FIG. 2: The distribution of ground state energies Eo 
{left panel) and crystal cell size {right panel). The 
colors indicate the nature of the FLS. 


decrease in packing density. 

All but one of the racemic groundstates, for which 
both enantiomers are equally favoured, consist of actually 
racemic crystals with equal numbers of each enantiomer. 
This outcome is not imposed by the Hamiltonian, as evi¬ 
denced by the one exception, the racemic FLS 12r, whose 
groundstate contains only a single enantiomer, resulting 
in a crystal structure identical to that of the chiral option 
12c. This crystal is depicted in Figure[^a). In the case of 
chiral organic molecules, there is considerable interest in 
understanding the conditions under which crystallization 
might select a single enantiomer from a racemic mixture. 
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Starting with Pasteur’s famous separation of the enan¬ 
tiomers of tartaric acid m, enantiomer-selective crystal¬ 
lization represents an important route to obtaining single 
enantiomer populations of molecules. Empirically, chiral 
organic molecules are found to crystallize preferentially 
into the pure enantiomer in 10% of cases [20]. Our one 
example, the FLS 12r, is insufficient to establish any gen¬ 
eral rule. That said, we note the groundstate energy of 
the chiral crystal 12c, Eq = —3/4, is only case of a chiral 
crystal with an energy lower than —1/2. 

Our FLS model clearly demonstrates that the racemic 
crystal is generally more stable than its chiral counter¬ 
part. This stability is highlighted by the fact that, of the 
three choices of FLS which achieve perfect packing, i.e. 
Eo = —1, two are racemic; 74r and 75r (see Figure [^b) 
for the structure of the latter case). These perfect pack¬ 
ings, along with the general trend for enantiomers to pack 
in equal proportion point to a general packing benefit 
that is afforded low symmetry FLS’s when they can be 
paired with their enantiomer. The resulting dimer has 
a higher symmetry (gaining, at the least, a centre of in¬ 
version) and this allows for a greater number of crystal 
types to be considered which, in turn, would be expected 
to lead to a lower groundstate energy. 


B. Crystalline cell sizes 

The unit cell size Z of the groundstate crystal struc¬ 
ture provides some measure of the structural complexity 
of the crystal since the larger Z, the number of lattice 
sites within the unit cell, the greater the information re¬ 
quired to specify the structure. In Figure we plot the 
distribution of unit cell sizes, revealing a large variation, 
ranging from 1 to 128, with peaks at powers of 2. A se¬ 
lection of crystalline ground state structures is presented 
in Figure 

The size of the largest unit cells observed in the FLS 
model is of particular interest. On a practical point, a 
unit cell of 128 sites highlights the diffiiculty of ensuring 
that the simulation volume is appropriately shaped to ac¬ 
commodate the crystal. The occurrence of crystals with 
large unit cells is also a topic of considerable physical in¬ 
terest. Pauling, in 1955, solved the structure of an inter- 
metallic alloys, NaCd 2 , to find that the unit cell was huge, 
containing 1152 atoms [21]. Since this pioneering result, 
intermetallic alloys have provided a growing number of 
what have been labelled giant unit cell structures [22] . 
with hundreds of cases of cells whose size exceeds 10^ 
atoms. The study of such giant unit cell structures has 
been identified [22] as ’’one of the last white spots on the 
map of metal physics”. At issue are the twin problems of 
the description of their organization and of the identifi¬ 
cation of the mechanism by which they are formed from 
the liquid state. 

We find two chiral FLS’s, 28c and 67c, which exhibit 
groundstates with a cell size Z = 128. The groundstate 
of the the FLS 67c is shown in Figure [^e). Despite 


this apparent structural complexity, we find both crystals 
nucleate and grow spontaneously in Monte-Carlo simu¬ 
lated annealing of sufficiently large systems, confirming 
that they are indeed kinetically accessible from the super¬ 
cooled liquid. A preliminary examination of the crystal 
structures of the 28c and 67c groundstates indicate the 
presence of asymmetric unit cells, 8x4x4 and 8 x 8 x 2 re¬ 
spectively, characterised by repeated motifs of clusters of 
FLS’s. Unravelling the energetics and kinetics of assem¬ 
bly of such hierarchical structures remains an interesting 
challenge. 


IV. THERMAL PROPERTIES 

We now turn to the study of the FLS model at finite 
temperature in the canonical ensemble. We first present 
the properties of the equilibrium liquids and the equi¬ 
librium freezing transitions. We then move to the non¬ 
equilibrium properties, by investigating the dynamics in 
the supercooled liquid and the low-temperature fate of 
an annealed system. 


A. Equilibrium liquids and freezing transitions 

Top panels of Figure [^present the temperature depen¬ 
dence of the average energy (ie., minus the FLS density) 
in canonical ensemble simulations of the model, for vari¬ 
ous choices of the FLS. A striking feature of these curves 
is the strongly first-order nature of the freezing transi¬ 
tions, which is generic of the model. In an hysteresis 
cycle, we observe both strongly supercooled liquid and 
overheated crystal. In order to locate the equilibrium 
freezing point, we employed a Maxwell construction, as 
in the 2D model [8]. The resulting temperature is indi¬ 
cated as a thick red vertical line in Figure 

1. Order and its Entropy Cost in the Liquids 

At infinite temperature, the liquid has a finite energy 
which we denote Fqo, due to random appearances of the 
local structure in a completely disordered spin system. 
These infinite temperature FLS densities do not exceed 
0.012. We can evaluate the ordering of the liquid by 
relating its FLS density to the asymptotic regimes T = 0 
and T = oc. We introduce the order parameter 0, defined 
as [7] 

0(T) = {E{T) - E^)/{Eo - F;oo) (1) 

In Figure [^ we present the distribution of the maxi¬ 
mum value of (j) attained in the equilibrium liquid {i.e. 
the value at the freezing point Tf). We find that it is 
very small in most cases, with an average of {(j){Tf)) = 
0.11(0.1), indicating negligible ordering in the liquid. It 
is less than 0.1 in 71 cases out of 115, and more than 0.3 
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FIG. 3: Top row: The energy versus temperature curves for a selection of FLS’s. Monte-Carlo simulations of an 
hysteresis cycle are shown in blue and green dots (respectively cooling and heating), with the dashed line indicating 
freezing and melting transitions. The thick red curve corresponds to a thermodynamic calculation of the equilibrium 
curve by a Maxwell construction. Bottom row: the energy autocorrelation time, in log-scale. 


in only five cases. This is in sharp contrast with the 2D 
model, for which equilibrium liquids with > 0.8 

were observed. 



Order parameter {E-E^ )/{E,-E^ ) 


FIG. 4: Statistics for the amount of ordering. Dark 
gray: at the equilibrium crystallization transition, and 
light gray: at the limit of stability of the supercooled 
liquid. 


The energetic gain of increasing the density of favoured 
local structures competes, at finite temperature, with the 
entropy loss associated to this increased order. This en¬ 
tropy cost can be assessed by performing a high tempera¬ 
ture expansion of the energy-entropy relation, as already 
discussed in the 2D case [5] , 

S = So.-^{E^-Ef + 0{{E^-Ef) ( 2 ) 

where Soo = ln(2) (the maximum entropy of a binary 
spins system), and A is the entropy cost parameter quan¬ 
tifying the decrease in entropy with decreasing energy. 
The larger the value of A, the more rapidly the liquid 


loses entropy as the density of FLS’s increases and the 
higher the temperature at which the liquid becomes ther¬ 
modynamically unstable with respect to the solid. We 
established in [8] that this high temperature expansion 
is actually a cluster expansion, where the coefficient A 
is related to the statistics of pairwise overlaps between 
neighbouring FLS’s. We calculated the values of A for 
the 3D FLS model and find them to be of the order of 
150, considerably larger than the analogous values for the 
2D model (which are of the order of 5). Thus the entropy 
cost for a specific FLS in 3D is considerably higher than 
that found in 2D. This result is due to the much larger 
number of local structural possibilities in 3D as compared 
to 2D and, hence, the greater entropy loss when a single 
FLS is selected. 

This high entropy cost places a serious limit on the den¬ 
sity of FLS’s the liquid state can physically accumulate. 
Kauzmann [23] noted that extrapolations of the liquid 
entropy to low temperatures suggested that the entropy 
would vanish at some non-zero temperature. Equation 
implies an analogous disappearance of the liquid entropy 
(and, hence, the ’end’ of any liquid state) at an energy 
Ek, related to A by 

Ek = E^- (3) 

Based on this , admittedly simplistic, estimate, the values 
of Ek range between —0.16 and —0.02, underlining the 
very limited amount of structure the 3D FLS liquids can 
accumulate. As discussed below, some of the supercooled 
liquids in the FSL model achieve energies much lower 
than —0.16 indicating the limits of this high-temperature 
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estimate, which neglects cooperative effects of more than 
2 FLS. 


2. Equilibrium Freezing Temperatures 

Given the low densities of the FLS in the equilibrium 
liquid, it is reasonable to approximate the liquid with its 
high temperature limit. Likewise, the crystal shows little 
decrease in order prior to melting, as shown in Figure 
Neglecting order in the liquid and loss of order in the 
crystal, we obtain an approximate formula for freezing 
temperatures: 


This expression depends solely on the ground state en¬ 
ergy, with no dependency on the liquid properties. In 
Figure we plot the freezing temperatures obtained by 
Maxwell construction against the ground state energies 
and compare them with the approximate expression of 
Equation!^ The data for the 2D model are also shown 
in the same graph. We find that Equation provides 
excellent agreement with the values of the freezing tem¬ 
perature for the ELS model on the ECC lattice. This 
success is in stark contrast to the 2D case where the the¬ 
oretical expression substantially overestimates Tf. 



EIG. 5: The freezing temperature Tf against the 
ground state energy for both the 3D (crosses) and 
the 2 D (circles) ELS model. The line corresponds to 
theoretical Equation 


B. Supercooling the Liquid 

Many of the ELS’s result in liquids that exhibit a 
marked supercooling before freezing. In this Section and 
the next, we employ a standard cooling rate of 5.10^ 
rejection-free MG steps per site and per unit of temper¬ 
ature, for system sizes of order 24^. 


1. General Charaeteristies of the Supereooled Liquid 

We find that the average ratio of the minimum tem¬ 
perature Tmin to which the liquid can be cooled to that 
of the freezing point to be (Tmin/Tf) = 0.79(0.09). This 
supercooling is significant and considerably larger than 
observed in the 2D model [ 8 ], apart from the one chi¬ 
ral case [9]. In molecular liquids, the glass transition Tg 
is typically associated with a deeper supercooling, i.e. 
Tg/Tf 0.66 [H]. 

Supercooling is associated with an increase in the den¬ 
sity of ELS above that observed at coexistence. In Eig- 
ure we report the distribution of (f calculated at the 
lowest value of the temperature to which each liquid can 
be supercooled. Here we find a handful of liquids reach¬ 
ing over half of the order of the groundstate crystal. The 
average is (0(T^m)) = 0.21(0.17), twice that at the equi¬ 
librium crystallization transition. 


2. Supereooling Scenarios 

In Eigurej^we plot the average energy E{T) and the 
energy autocorrelation time t(T) (measured in Monte- 
Garlo steps per site), estimated by the method of batch 
means [25] as a function of T. We have distinguished 
several different scenarios based on the behaviour of the 
temperature dependence of E and r in the supercooled 
regime. Examples of each scenario are presented. The 
properties of each scenario are described as follows. 

i) Non-Ordering Liquid (79 ELS’s). The most common 
scenario in this model is characterized by a liquid showing 
negligible cooperativity down to the freezing transition, 
as illustrated on Eigurej^. Even at the lowest tempera¬ 
ture Tmin at which the metastable liquid can be observed, 
the energy autocorrelation time remains small, only 3.5 
MG steps per site. 

ii) Ordering Liquid (22 ELS’s). The cases not included 
within the ’Non-Ordering’ category exhibit a significant 
accumulation of ELS. Eigurej^ is a good example of this 
scenario. While the equilibrium liquid exhibits almost no 
ordering {(j)(Tf) = 0.07), the supercooled liquid is char¬ 
acterized by a steady accumulation of ELS, with a high 
and approximately constant heat capacity Cy = dE/dT. 
In this regime, the dynamical time scale grows steadily 
up to r ~ 30 MG steps per site, indicating the existence 
of cooperative effects in the supercooled liquid. 

hi) Quasi-Critieal Liquid (14 ELS’s). In a number of 
cases, the relaxation time r is strongly non-Arrhenius 
and can climb up to 10 ^ MG steps per site (our limit of 
time resolution for r), exhibiting an apparent divergence 
at the end of the supercooled branch, with a correspond¬ 
ingly strong increase in heat capacity. This behaviour is 
illustrated on Eigurej^. The autocorrelation time is well 
fitted by a diverging power law r oc |T — with an 

exponent 7 ^ 1.3 (dashed line in Eigurej^). While this 
scenario resembles that of a glass transition, we note that 
this supercooled branch terminates with a discontinuous 
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transition in the energy in all cases in this model. 


C. Low Temperature States Following Quenching 

While we always observe a discontinuity in energy at 
the freezing transition when annealing the system, the 
resulting low-temperature state varies and the ground 
state crystal is not always obtained. We now discuss the 
different classes of T = 0 structure observed following a 
quenching via Monte Carlo dynamics. 

i) Standard Crystallization (99 FLS) The most com¬ 
mon scenario corresponds to the formation of a crystal 
structurally similar to the ground state computed from 
the enumerative method. There is some variability be¬ 
tween cases in this family. For 66 FLS’s, the observed 
state corresponds exactly to the computed ground state, 
provided that the boundary conditions were adapted to 
the crystal structure. In 9 FLS’s, we find the predicted 
groundstate crystal is formed but as polycrystalline, with 
domains separated by grain boundaries, as illustrated in 
Figure]^. In a rather large number of cases (24 FLS’s), 
the system cooperatively orders into a state with the 
same energy as the calculated ground state, but with 
apparent degrees of freedom remaining. We tentatively 
qualified these states as “plastic”, without intending to 
assert any specific connection with the more clearly de¬ 
fined plastic crystal phases of molecular systems [27] . 
This is reminiscent of the ground state structure of the 
antiferromagnetic Ising model on a triangular lattice, 
for which frustration effects lead to a highly degenerate 
ground state. Finally, we observe 3 FLS’s to crystallize 
into the predicted groundstate but via an intermediate 
mesophase whose stability remains unclear. An exam¬ 
ple of this multi-step crystallization is illustrated in Fig¬ 
ure p. Note that the mesophase is observed only on 
cooling of the liquid, and at temperatures lower than 
the computed equilibrium transition between liquid and 
crystal. 

We also identified two scenarios for which the low- 
temperature state was structurally and energetically dif¬ 
ferent from the predicted ground state. 

ii) Metastahle Polymorph Formation (14 FLS) In these 
cases, a recognizably crystalline structure was consis¬ 
tently obtained, albeit with an energy strictly higher 
than the ground state. This scenario is illustrated in 
Figure [^, for which a complex crystalline state with en¬ 
ergy —4/9 and a huge cell size Z > 288 was obtained, 
while a much simpler and stabler crystalline structure 
with energy —1/2 and cell size Z = 16 was available. 
In some cases the crystal state formed can vary between 
different cooling runs. This is clear evidence of polymor¬ 
phism. We do not know, in each case, whether the higher 
energy structure is, in fact, a true equilibrium phase at 
the freezing temperature, favoured over the lower energy 
structure by entropic effects, or a metastable crystal, se¬ 
lected by the kinetics of ordering. The demonstration 
that our model includes polymorphs is interesting given 



T E^=-2/5=-0A E^^ = -Q.3A9 

(a) Defected crystal formation in FLS 43r. 



T =-l/2=-0.5 =-4/9;^-0.44 

(b) Metastable polymorph formation in FLS 39 



^ 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Grouncl State Monte-Carlo 


T =-3/10 =-0.3 r;^i^ = -0.264 

(c) Amorphous solid formation in FLS 51c. 

FIG. 6: Three scenarios for avoided crystallization. In 
each case, we display: (i) the E{T) curve with a picture 
of the FLS in inset and (ii) the energy autocorrelation 
time r(T), as in Fig. (hi) the computed crystalline 
structure and (v) a low-temperature state from 
Monte-Carlo simulations; (iv) and (vi) the respective 
associated location of the favoured local structures 
(FLS sites in blue, non-FLS in white). 












































the prevalence of polymorphs in molecular crystals. In 
one well studied case [26], a molecule has been observed 
to crystallize into 7 distinct crystal structures. 

hi) Amorphous Solid Formation Finally, in a single 
case, the FLS 51c, the liquid freezes with an apparent 
discontinuity of the energy into an amorphous state with 
no recognizable trace of crystallinity (see Figurej^). This 
transition is accompanied by a dynamical arrest of the su¬ 
percooled liquid, in the quasi-critical scenario discussed 
in Section |IVB[ and the relaxation time scale becomes 
larger than our resolution of r = 10^ MC steps per site. 
On heating this state, negligible hysteresis is observed. It 
is not clear whether this puzzling behaviour corresponds 
to a conventional glass transition (for which the energy 
would have to be continuous) or some sort of polycrys¬ 
talline aggregation. What is clear is that the 51c case 
represents the interesting possibility of the formation of 
an amorphous solid stabilized entirely (by construction) 
by the same FLS as stabilizes the crystal. 


V. ON THE CORRELATION BETWEEN THE 
FLS STRUCTURE AND THE PROPERTIES OF 
THE LIQUID AND SOLIDS STATES 

In the previous sections, we have described the vari¬ 
ety of behaviours of our model, going through different 
aspects of its macroscopic properties: crystalline ground 
states, liquids and supercooled liquids, freezing transi¬ 
tions and low-temperature states. In this Section, we ex¬ 
amine a number of hypotheses associated with how the 
choice of FLS determines the collective properties of the 
system. A summary of our conclusions and the statistical 
data on which they are based is presented in Table [I| 

A. Structure of the Ground State 

We start by considering the influence of the symme¬ 
try properties of the FLS (the size S of the spatial point 
group of an FLS introduced in Secti on [ll| )on the struc¬ 
ture of the ground state. In Section m we introduced 
two observables to quantify the latter: the ground state 
energy Eq measures the quality of the packing, while the 
cell size Z measures its complexity. We disentangle the 
correlations between these quantities by restricting the 
statistical analysis to sets of FLS for which one parame¬ 
ter is fixed - for instance, we analyze the correlation be¬ 
tween S and Eo by considering the 22 FLS’s with crystal 
cell Z = 16. This is made possible by the existence of 
sharp peaks in the histograms of Figures Because Z 
and S vary over several orders of magnitude, we consider 
their logarithm when measuring correlations. 

The connection between symmetry and ground state 
energy is investigated in Table There appears to be 
no significant correlation between Eq and S: the degree 
of symmetry of the local structure does not influence its 
ground state energy. This is perhaps surprising, as one 


might expect awkward objects of low symmetry to have 
poor packing properties. We interpret the independence 
between Eq and S as the compensation of two effects: on 
the one hand, high-symmetry objects are more likely to 
fit well together; on the other hand, low-symmetry ob¬ 
jects have more rotational variants, thus offering a wider 
pool of combinations to select the ground state from. 

On the other hand, S and Z are significantly corre¬ 
lated, with a correlation ratio r = —0.64, as reported 
in Table ib. Thus the cell size increases with decreas¬ 
ing symmetry. However, this trend is loose and admits 
exceptions: for instance, highly symmetric FLS 7r, for 
which 5 = 6, has a large unit cell Z = 64 depicted 
in Figure [^c). On the other hand, low-symmetry FLS 
32c (with 5 = 1) crystallizes in a simple structure with 
Z = 7, as depicted in Figure [^d). 

B. The special role of the chiral symmetry 

As already hinted in Section |III A[ the ground state 
energies of chiral FLS are generally higher than those of 
achiral structures. This claim is properly confirmed in 
Table U, by comparing two groups of FLS with identical 
degree of symmetry: on the one hand the 19 achiral FLS’s 
with a plane reflection symmetry only, and on the other 
hand the 10 chiral FLS’s with only a rotational symmetry 
of order 2. We observe a statistically significant differ¬ 
ence in the average ground state energy, while the differ¬ 
ence in crystal cell size is insignificant. This shows that 
plane reflection symmetry of the structure plays a spe¬ 
cial role in determining the ground state energy: chiral 
local structures pack less densely than achiral structures. 
This special role of the chiral symmetry in such discrete 
packing problems is, to our knowledge, an original result 
of this study. 


C. Crystallization or not? 

We next turn to the question of the influence of struc¬ 
tural properties on the low temperature fate of the sys¬ 
tem. We consider on the one hand the 90 “crystallizers” 
which find their exact ground state, and on the other the 
16 FLS’s which do not find their ground state structure, 
and form either a glassy state or a metastable polymorph. 
We exclude defected crystals from this analysis, as they 
lie in-between. The dominant contribution is the sym¬ 
metry factor: low symmetry strongly conditions towards 
avoided crystallization, as demonstrated in Table p. 

Because low symmetry structures also tend to form 
large crystal cells (Table crystal formation also 

strongly correlates with crystal cell size. However, we 
can refine this observation by considering the set of FLS 
with 5 = 1. In Table [l^, we show that the value of Z 
does not significantly affect the probability to fail crys¬ 
tallization within this population. To conclude, we have 
found that low symmetry FLS tend to crystallize less 
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TABLE I: Summary of our correlation hypotheses and associated conclusions. Their validity is assessed by the 
p-value, corresponding to the likelihood of observing comparable correlations in the case of independent variables, as 
provided by either Pearson’s analysis* or Student’s t-test’*'. Numbers in parentheses are standard deviations. 


Question 


Conclusion 


Method employed and results 


(a) Does the symmetry of the 
ELS affect the ground state 

energy? 

(b) Does the symmetry of the 
ELS affect the crystal 

structure? 

(c) Does chirality of the ELS 
inffnence the ground state 

energy? 

(d) Does the symmetry of the 
ELS affect the likelihood to 

form the ground state crystal? 

(e) Does the crystal cell size 
further affect the likelihood to 
fail crystallization? 

(f) Is there a link between the 
supercooling scenario and the 
gronndstate properties? 


Compare Eo and log<S at fixed Z 


Data set 

points 

Pearson r 

p-value* 

Z = 16 

22 

0.015 

0.95 


Compare log Z and log S at fixed Eo • 


Data set 

# points 

Pearson r 

p-value* 

Eo = -1/2 

35 

-0.64 

10“® 


Compare Eo for chiral and achiral populations. 


Data set 

# points 

(Eo) 

p-value^ 

Achiral, 5 = 2 

19 

-0.55 (0.10) 

0.003 

Chiral, 5 = 2 

10 

-0.41 (0.14) 



Correlate low temperature state to symmetry 

Oa.tfl. spt zkk nnints ( Ti-^T^alv 


Data set 

# points 

(5) 

p-value^ 

Crystallizing FLS 

90 

3.0 (5.3) 

0.007 

Avoided crystallization 

16 

1.25 (0.4) 



Correlate low-T state to Z at fixed symmetry S. 


Data set 

# points 

{E) 

p-value^ 

Crystallizing FLS, 5=1 

34 

31 (27) 

0.2 

Avoided cryst., 5 = 1 

12 

38 (22) 



Correlate the supercooling scenario to Z and Eo. 


Non-ordering 
{Eo) -0.56 (0.17) 
(Z) 18 (18) 


Ordering 
-0.42 (0.12) 
32 (21) 


Quasi-critical p-val.* 

-0.36 (0.08) 10"® 
50 (32) 10"^ 


No, they are nncorrelated. 


Yes, Z increases as S decreases. 


Yes, chiral ELS have a higher 
Eo than achiral ones. 

Yes, low symmetry ELS’s are 
more likely to avoid the 
gronndstate. 

No, Z does not inffnence 
crystallization at fixed S. 

Yes, the supercooling scenario 
is linked to Z and Eo (see 
discussion). 


than high symmetry ones, but that the crystal cell size 
- and hence its “complexity” - is no further impediment 
to crystallization. We note in particular that our three 
“giant unit cell” crystals with cell size Z > 80 straightfor¬ 
wardly freeze into a perfect crystalline structure. On the 
other hand a low ground state energy favours crystalliza¬ 
tion. The latter observation can be intuitively explained 
by noting that a crystalline structure that is a deep min¬ 
imum of energy - with low Eo - is less likely to have 
viable competing polymorphs than a high energy ground 
state. 

Only a single ELS leads to an amorphous solid, as 
discussed in Section |IV C[ and we cannot make general 
statements on this basis. We note however that it is chi¬ 
ral, with minimal symmetry 5 = 1, has a high ground 
state energy Eo = —3/10, and a large crystal cell Z = 80. 
Variants or extensions of the present model could provide 
a larger population of similar amorphous systems, and 
give a statistical basis to draw more general trends. 


ordering, and divergent time scale (“quasi-critical”). In 
Table [J we correlate this classification with the ground 
state energies and crystal cell sizes of all ELS’s. 

We observe a clear trend: high order in supercooled liq¬ 
uids correlates with high ground state energy and large 
crystal cells. The first observation is not surprising. In¬ 
deed, as we have seen in Section [IV A[ high ground state 
energy correspond to low freezing temperatures, and thus 
to the possibility for the liquid to overcome high entropic 
barriers to ordering. On the other hand, the fact that 
crystal structure correlates with liquid ordering is rather 
unexpected, as in a simple nucleation picture the liquid 
does not “know” about the crystal structure. This is not 
a symmetry effect: restricting to the 26 chiral structures 
with 5 = I, we observe the same effect (although not 
resolved statistically). 


VI. DISCUSSION 


D. Supercooling and ground state structure 


The last point of our statistical analysis concerns the 
physical properties of the supercooled liquid and its con¬ 
nection with structural properties. In Section |IVB[ we 
qualitatively distinguished three different behaviours in 
the supercooled liquid: no trace of ordering, significant 


What are the macroscopic consequences of local struc¬ 
ture in liquids? In this article, we addressed this ques¬ 
tion in the case when the favoured local structures found 
in the liquid and in the solid are the same - that is, 
when the ground state corresponds to a dense accumu¬ 
lation of a single ELS. We introduced a simple three- 
dimensional lattice model in which the ELS is explicitly 
specified. This Eavoured Local Structures model on a 3D 
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face-centered cubic lattice provides a discrete set of 115 
such structures to choose from, and we reported here on 
the comprehensive study of these systems. This large li¬ 
brary of structures is a unique feature of this work and 
allowed us to employ an original statistical approach to 
correlate observable features of the model and unfold un¬ 
derlying physical laws. Given the similarity between dif¬ 
ferent instances of the model - only the geometry of the 
FLS varies - the variety of observed behaviours is note¬ 
worthy. 

We focused on the structural properties of the ground 
states - studying the relation between the symmetry of 
an FLS and its packing properties - and on the proper¬ 
ties of the supercooled liquid - addressing the questions 
of how much order can be accumulated prior to freezing, 
and what solid phase will be formed at the freezing tran¬ 
sition. Our main results can be summarized as follows: 

1. All ground states are crystalline, with high variability 
in both the ground state energy —1 < Eq < —1/4 and 
unit cell size 1 < ^ < 128. This structural complex¬ 
ity does not appear to present any systematic kinetic 
obstacle to crystallization. 

2. We found no evidence that low symmetry FLS’s re¬ 
sult in ’poor’ {i.e. high energy) crystals. Instead, we 
found that low symmetry FLS’s tend to produce more 
complex crystals as characterised by larger unit cells. 
We interpreted this as resulting from compensating ef¬ 
fects: low symmetry objects offer more orientational 
freedom than high symmetry ones, but the latter gen¬ 
erally fit better together. We observed an exception 
to this rule in the case of chiral structures, which pack 
poorly, compared to achiral structures with equal de¬ 
gree of symmetry. 

3. We observed only a single case of enantiomeric sep¬ 
aration of a racemic mixture by crystallization. This 
exception corresponding to a chiral FLS that managed 
a particularly low energy groundstate. 

4. The entropic cost of ordering in the liquid is so high 
in our model that we observe little organization prior 
to equilibrium crystallization. We regard this a conse¬ 
quence of the very strict definition of local structures 
in our model. Considering incomplete or multiple FLS 
could allow for more complex liquid behaviour. 

5. Freezing transitions are typically strongly first-order 
and most liquids can be readily supercooled. Ap¬ 
proaching the freezing transition, some liquids ex¬ 
hibit significant ordering and even apparently diverg¬ 
ing relaxation times, while others remain unstruc¬ 


tured. High ground state energy depresses the transi¬ 
tion point towards low temperatures, allowing higher 
order in the liquid. More surprising, we found that sys¬ 
tems whose crystalline structure is complex (i.e. large 
Z) exhibit stronger ordering in the supercooled liquid, 
showing a significant correlation between the ground 
state structure and the amount of order in the amor¬ 
phous state. 

6. Finally, we found that low symmetry and high ground 
state energy FLS tend to avoid crystallization by ei¬ 
ther falling into a metastable crystalline structure, or 
into an amorphous state. On the other hand the com¬ 
plexity of the crystalline structure does not affect the 
likelihood to crystallize. In particular, our giant unit 
cell crystals are readily formed in Monte-Carlo simu¬ 
lations of large systems. 

Perhaps the most telling conclusion from this study 
comes from a negative result. Despite our comprehen¬ 
sive inclusion of all possible local structures, we never 
found a local structure so geometrically ’awkward’ that 
it lacked a crystalline groundstate: our enumerative pro¬ 
cedure to find crystalline structures always provided the 
lower bound on observed energies. Furthermore, every 
choice of FLS (with a single exception), exhibited a rel¬ 
atively rapid crystallization, if not always to the ground 
state. While the choice of FLS can influence the size 
of the unit cell in the crystal and the degree of order¬ 
ing in the supercooled liquid, we conclude that the exis¬ 
tence and kinetic accessibility of the crystalline phase are 
essentially independent of the nature of the local struc¬ 
ture. As for the liquid state, it is clearly the multiplic¬ 
ity of the FLS, not its geometrical structure, that exerts 
the most important influence with regards the accumula¬ 
tion of structure in the liquid and any associated slowing 
down. This conclusion challenges the long lived tradition 
of assigning significance to the presence of specific local 
structures in low temperature liquids. 
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Appendix A: Enumerative method for determining 
ground states 

Our method for determining ground state structures 
and energies makes use of a correspondence between crys¬ 
talline base vectors on a lattice and number-theoretical 
objects called Hermite matrices. A similar approach was 
introduced by Hart and Forcade m as a tool to study 
derivative structures for metallic alloys. The idea is to 
list all possible sets of crystalline base vectors (superlat¬ 
tice) on the FCC lattice. There are 17602 distinct super¬ 
lattices of the FCC lattice with up to Z = 80 sites per 
crystalline cell (having eliminated redundancies due to 
reparametrization and lattice symmetries). We use the 
correspondence between superlattices and integer matri¬ 
ces in the Hermite normal form, completed by a reduction 
using the FCC symmetries, to obtain this set of struc¬ 
tures. For each such superlattice and every FLS, we find 
the best crystalline spin configuration associated to this 
superlattice by simulated annealing of the spin variables 
of a single crystal cell, with appropriate boundary condi¬ 
tions. The small number of spins involved in this makes 
us confident that the best spin structure is found. When¬ 
ever necessary, we extend this procedure up to 128 sites 
per cell. This systematic approach is complemented with 
a more standard Monte-Carlo simulated annealing using 
a rejection-free algorithm El of a large system of up to 
24^ sites, which always gives a resulting minimal energy 
higher or equal to that obtained in the enumerative pro¬ 
cedure. While these approaches can in principle fail and 
fall in metastable minima, we are confident that most (if 
not all) of the ground states are indeed the global ones. 


